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I.  Introduction 

This  report  describes  an  empirical  analysis  of  a  motion 
model  that  has  Been  used  to  generate  random  submarine  tracks  for  an 
ASW  tactical  decision  aid.  The  model  describes  a  submarine’s 
motion  as  a  series  of  transitions  between  the  square  cells  of  a 
grid  that  covers  a  defined  operating  region.  A  3  X  3  transition 
matrix  is  associated  with  each  cell  of  the  grid  which  determines 
the  submarine's  transitions  from  a  cell.  The  set  of  transition 
matrices  define  a  Markov  process.  Despite  its  discrete  nature, 
this  Markov  track  generating  process  has  been  called  a  diffusion 
process  in  antisubmarine  warfare  tactical  decision  aid  literature. 
The  transition  matrices  are  determined  by  tracks  generated  by  an 
auxiliary  stochastic  process  that  is  presumed  to  be  of  higher 
fidelity  but  more  costly  to  implement  than  the  Markov  process. 

The  auxiliary  track  generating  process  that  was  used  in  the 
analysis  is  a  random  tour  process  that  was  developed  to  generate 
operationally  realizable  submarine  tracks.  The  process,  which  is 
described  in  detail  in  Reference  1,  generates  random  tracks  that 
have  specified  end  points  and  a  specified  length  or,  for  a 
submarine  with  a  constant  track  speed,  a  specified  transit  time. 
The  auxiliary  process  is  implemented  by  the  program  that  is  listed 
in  Appendix  1. 

The  auxiliary  tracks  are  the  basis  for  the  generation  of  two 
types  of  diffusion  transition  matrices  by  a  procedure  described  in 
Section  II.  The  first  type,  called  static  transition  matrices,  are 
the  transition  matrices  that  define  the  subject  Markov  motion 
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model.  A  static  transition  matrix  is  associated  with  each  cell  of 
the  operating  region  grid.  The  second  type,  called  dynamic 
transition  matrices,  define  an  alternate  Markov  motion  model. 

This  model  differs  from  the  subject  model  in  that  the  transition 
matrix  associated  with  a  cell  depends  on  the  number  of  transitions 
that  have  taken  place  in  the  motion,  that  is,  it  is  time  dependent. 

The  auxiliary  tracks  are  also  the  basis  for  the  generation  of 
a  set  of  maps  that  describe  a  submarine's  position  at  a  sequence  of 
equally  spaced  times  that  are  determined  by  a  fixed  time  step.  The 
number  of  times  a  cell  is  occupied  at  the  end  of  a  time  step 
divided  by  the  number  of  tracks  is  an  estimate  of  the  probability 
that  the  submarine  was  within  the  cell  at  the  time  step.  For  this 
reason,  such  maps  are  called  probability  maps  in  the  following 
discussion. 

Tracks  generated  by  the  subject  static  diffusion  Markov 
process  and  those  generated  by  the  alternate  dynamic  diffusion 
Markov  process  also  can  be  used  to  generate  probability  maps.  Some 
differences  between  the  three  track  generating  processes  is  evident 
by  a  comparison  of  the  random  tour  probability  maps,  the  dynamic 
diffusion  probability  maps  and  the  static  diffusion  probability 
maps  that  are  shown  on  the  following  pages. 
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li.  The  Transition  Matrix  Generation  Procedure 

The  transition  matrices  that  are  used  in  the  analysis  were 
determined  by  the  program  listed  in  Appendix  1.  In  the  program,  a 
random  tour  track  is  generated  by  a  submarine  moving  with  constant 
speed  between  two  points  located  in  a  plane.  The  points  lie  on  the 
x-axis  of  a  rectangular  coordinate  system  whose  origin  is  midway 
between  the  points.  Lines  parallel  to  the  x-axis  and  lines 
parallel  to  the  y-axis  define  a  grid  of  square  cells  in  the  plane, 
(For  examples  of  the  program  geometry,  see  Section  V.)  In  the 
following  discussion  of  the  transition  matrix  determination 
procedure,  grid  cells  are  identified  by  the  rectangular  coordinates 
of  their  centers. 

During  its  motion,  the  cell  occupied  by  the  submarine  is 
determined  at  a  sequence  of  times  separated  by  a  fixed  time  step. 
From  this  determination,  transition  matrices  are  constructed  as 
follows:  First,  the  cell  that  is  occupied  is  determined  for  each 
time  step.  For  I  =  1  to  NS  (the  number  of  time  steps) ,  if  cell 
(LO,  MO)  is  occupied  by  the  submarine  at  the  end  of  the  I-l®^  time 
step,  at  the  end  of  I**’  time  step,  the  submarine  will  be  in  a  cell 
(LO  +  N,  MO  +  K)  where  N  =  -1,  0,  or  1  and  K  =  -1,  0,  or  1.  In 
this  case,  1  is  added  to  the  (N,K)  element  of  a  3  X  3  transition 
matrix  associated  with  cell  (L0,M0)  for  time  I-l.  The  elements  of 
this  matrix  are  stored  in  an  array  MT(N,  K,  LO,  MO,  I-l).  In 
addition,  1  is  added  to  the  (N,K)  element  of  a  second  3X3 
transition  matrix  associated  with  cell  (L0,M0).  And,  the  elements 
of  this  matrix  are  stored  in  an  array  MTS(N,  K,  LO,  MO).  By 
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repeating  this  procedure  for  each  random  tour  track,  and  dividing 
the  resulting  matrix  elements  by  the  number  of  tracks,  the  elements 
of,  the  3X3  matrix  in  the  array  MT  become  estimates  of  the 
elements  of  the  dynamic  transition  matrix,  and  the  elements  of  the 
matrix  of  the  3X3  matrix  in  the  array  MTS  become  estimates  of 
the  elements  of  the  static  transition  matrix. 

The  transition  matrices  and  the  probability  maps  used  in  the 
analysis  were  generated  from  10,000  random  tour  tracks.  Since  the 
time  step  duration  is  equal  to  the  cell  side  length  divided  by  the 
track  generating  submarine's  speed;  during  a  time  step,  the 
submarine  will  either  transition  to  an  adjacent  cell  or  remain  in 
its  current  cell.  This  results  in  a  transition  matrix  of  nine 
elements,  each  element  giving  the  probability  of  a  transition. 
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III.  Some  Compef Isons  of  the  Track  Generating  Processes 

The  dynamic  and;  static  transition  matrices  of  the  two  discrete 
Markov  track  generating  processes  generate  tracks  that  are  random 
but  differ  in  fine  structure  from  those  generated  by  the  random 
tour  process.  This  is  not  surprising  considering  the  loss  of 
information  that  occurs  in  their  generation.  Since  the  Markov 
motion  is  from  a  cell  center  to  the  cell  center  or  a  neighboring 
cell  center,  only  nine  motions  are  possible.  In  particular,  for  a 
submarine  that  moves  with  constant  track  speed,  this  results  in  a 
loss  of  position  and  velocity  information.  Relative  to  the 
auxiliary  track  generating  process,  the  diffusion  tracks  no  longer 
satisfy  the  constraints  of  constant  speed,  specified  length  and 
specified  final  end  point. 

The  fidelity  of  the  diffusion  tracks  could  be  improved  in  two 
ways;  First,  by  reducing  the  duration  of  the  time  step  and 
consequently  the  size  of  the  cells,  less  positional  information 
would  be  lost.  Second,  by  generating  a  velocity  (direction  of 
motion)  distribution  for  each  cell  and  each  time  step,  transition 
matrices  could  be  constructed  that  would  be  more  descriptive  of  the 
random  tour  motion  and  less  directional  information  would  be  lost. 
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IV.  Conclusions 

The^  primary  conclusion  from  the  analysis  is:  For  discrete 
Markov  processes  whose  state  space  if  current  location  and  time, 
there  are  tactically  reasonable  submarine  motion  scenarios  that  can 
not  be  closely  approximated  by  tracks  generated  with  static 
transition  matrices  but  can  be  closely  approximated  with  dynamic 
transition  matrices.  There  are,  however,  other  tactically 
reasonable  ^submarine  motion  scenarios  that  can  not  be  closely 
approximated  with  either  procedure.  In  what  is  possibly  the 
simplest  example  of  this,  consider  a  submarine  that  moves  among 
three  adjacent  cells  labeled  1,  2  and  3.  Suppose  the  submarine 
chooses  either  path  (1,2,3)  or  path  (3,2,1),  each  with  probability 
.5.  Given  the  submarine  is  in  cell  2,  the  probability  that  the 
submarine  will  transition  to  cell  3  is  dependent  on  its  previous 
location  as  well  as  on  its  current  location.  Consequently,  a 
Markov  process  whose  state  space  is  only  current  position  and  time 
can  not  describe  its  motion. 

For  the  random  tour  scenarios  that  were  investigated,  dynamic 
diffusion  tracks  were  found  to  rather  closely  emulate  the  random 
tour  tracks  from  which  the  dynamic  transition  matrices  were 
derived.  In  addition,  and  in  contrast  to  the  static  diffusion 
tracks,  the  dynamic  diffusion  tracks  generated  probability  maps 
that  duplicated  the  random  tour  probability  maps.  Although  the 
static  tracks  were  generally  similar  in  appearance  to  the  dynamic 
diffusion  tracks,  they  did  not  uniformly  terminate  at  the  random 
tour  track  final  end  point.  However,  for  the  chosen  scenario. 
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approximately  half  of  the  tracks  did  dp  so  as  is  evident  from  the 
static  diffusion  probability  maps  that  are  shown  in  Section  V. 

The  divergence  between  the  dynamic  diffusion  probability  maps 
and  the  static  diffusion  probability  maps  is  even  more  evident  when 
the  initial  point  and  final  point  of  a  track  coincide.  The 
explanation  for  this  is  the  variation  of  the  dynamic  transition 
matrices  with  time.  That  is,  the  same  cell  must  perform  a 
"submarine  expansion"  function  early  in  the  motion  and  a  "submarine 
contraction"  function  towards  the  end.  But,  when  the  starting  and 
ending  cells  are  distant  from  each  other,  the  function  of  any  cell 
does  not  change  significantly  with  time.  That  is,  cells  near  the 
start  cell  are  for  always  for  "expansion"  and  those  near  the  end 
cell  are  always  for  "contraction." 

The  random  tour  process  on  which  the  empirical  analysis  is 
based  was  chosen  to  amplify  the  divergence  between  the  random  tour 
tracks,  the  dynamic  diffusion  tracks  and  the  static  diffusion 
tracks.  The  significance  of  this  divergence  for  a  tactical 
decision  aid  incorporating  a  motion  model  based  on  static 
transition  matrices  will  depend  on  the  desired  applications  of  the 
decision  aid.  For  example,  it  might  be  less  significant  if  the 
location  of  the  final  end  point  of  a  submarine  track  were  chosen 
randomly  or  if  computer  memory  or  processing  speed  limited  modeling 
choices. 
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The  figures  in  Section  V  that  follows  show  probability  maps 
and  tracks  generated  by  the  random  tour  process ,  the  dynamic  Markov 
process  transition  matrices,  and  the  static  Markov  process 
transition  matrices.  These  figures  illustrate  the  relative 
effectiveness  of  the  dynamic  and  static  Markov  motion  model  in 
emulating  the  random  tour  motion  model.  The  diffusion  tracks  that 
are  shown  were  generated  with  17  diffusions.  This  results  in 
tracks  that  are  approximately  equal  in  length  to  the  length  the 
random  tour  tracks. 
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V.  The  Analysis  Progreua  Output 

The  figures,  tables  and  maps  that  follow  are  based  on  data 
generated  by  the  analysis  program  that  is  listed  in  Appendix  1. 
They  provide  the  basis  for  the  empirical  analysis. 
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Figure  l.  A  track  generated  by  a  random  tour  process  used  in  the 
analysis.  The  larger  circles  on  the  left  and  the  right  are  the 
initial  and  final  point  of  the  track.  The  smaller  circles  are 
track  positions  at  step  times. 
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Figure  2.  A  superposition  of  10,000  random  tour  tracks.  These 
tracks  determined  the  Markovian  transition  matrices  that 
generated  the  diffusion  probability  maps  that  follow. 
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RANTRACK. BAS  random -tour  probability  map 
number  of  time  steps  =  15 
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Table  i.  A  dynamic  diffusion  transition  matrix  that  was  used  to 
generate  the  dynamic  diffusion  map  for  1  transition. 


I^TRACk.BAS  dynamic  transition  matrix 
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Figure  3  A  track  generated  by  a  dynamic  process.  The  dynamic 
transition  matrices  were  determined  by  the  random  tour  tracks 
shown  in  Figure  2.  The  number  of  diffusions  is  17.  The  circles 
are  defined  in  Figure  1. 
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Figure  4  A  superposition  of  10,000  dynamic  diffusion  tracks. 
Their  dynamic  transition  matrices  were  determined  by  the  random 
tour  tracks  shown  in  Figure  2.  The  number  of  diffusions  is  17. 
The  circles  are  defined  in  Figure  1. 
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Table  il.  A  static  diffusion  transition  matrix  that  was  used  to 
generate,  the  static  diffusion  maps. 
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Figure  5  A  track  generated  by  a  static  diffusion  process.  The 
static  transition  matrices  were  determined  by  the  random  tour 
tracks  shown  in  Figure  2.  The  number  of  diffusions  is  17.  The 
circles  are  defined  in  Figure  1. 
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Figure  6  A  superposition  of  10,000  static  diffusion  tracks. 
Their  static  transition  matrices  were  determined  by  the  random 
tour  tracks  shown  in  Figure  2.  The  number  of  diffusions  is  17. 
The  circles  are  defined  in  Figure  1. 
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>MNTRACI(;'BAS  statfC'dlffusion  proi^bility  map 
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yi.  The  Analysis  Progreun 

The  analysis  program  requires  the  following  inputs:  the  speed 
of  the  submarine  in  knots;  the  distance  between  the  end  points  in 
nautical  miles;  the  track  length  in  nautical  miles;  delta,  the 
average  distance  between  course  changes  (leg  length)  in  nautical 
miles.;  the  leg  length  distribution  parameter,  an  integer  that 
defines  the  gamma  distribution  that  determines  leg  length  ( Input ing 
1  gives  an  exponential  distribution.)  and  the  number  of  tracks  to 
be  generated.  If  an  auxiliary  pseudorandom  number  generator  is  not 
chosen,  the  option  to  supply  a  seed  for  the  QuickBASIC  random 
number  generator  is  given.  Next,  the  option  to  either  view  single 
tracks  and  terminate  the  program  or  to  input  the  number  of  tracks 
to  be  generated  to  produce  data  of  the  kind  presented  above. 

The  program  is  written  for  a  computer  monitor  with  screen  mode 
12  capability.  To  change  this  requirement,  the  statement  SCREEN  12 
on  Line  4780  must  be  changed.  For  example,  for  a  monitor  with  CGA 
capability,  it  could  be  replaced  by  SCREEN  2. 

To  run  the  program  under  QuickBASIC,  the  /ah  switch  must  be 
used  to  duplicate  the  program  results  in  this  report.  This  switch 
allows  arrays  of  size  greater  than  64K.  In  addition,  it  may  be 
necessary  to  have  FILES  =  8  and  BUFFERS  =2  in  the  CONFIG.SYS 
file  and  to  remove  TSR  file  activation  statements  from  the 
AUTOEXEC.BAT  file. 

The  program  code  listed  in  Appendix  1  provides  for  the 
addition  of  an  auxiliary  random  number  generator  as  an  alternative 
to  the  QuickBASIC  generator  through  Lines  200,  210  and  220.  These 
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lines  allow  auxiliary  generator  setup  code  in  a  subroutine  starting 
at  Line  7720  and  generator  code  in  a  subroutine  starting  at  Line 
7810.  This  line  number  was  chosen  to  accommodate  code  that  is 
based  on  a  Generalized  Feedback  Shift  Register  (GFSR)  pseudorandom 
number  generator  that  is  listed  in  Reference  3.  If  an  auxiliary 
generator  will  not  be  added,  Lines  200,  210  and  220  can  be  removed. 
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Appendix  1.  The  Analysis  Progriua  Listing 


10  N$  S' "RANfRACK.BAS":  REN  a  program  to  analyze  some  random  track  generation  procedures 
20  CIS  :  OEFINf  I-M. 

30  PI  s  4  *  ATN(I);  NTR  =  0;  R$  =  '.'standard'';  HaG3  =  0;  FLAG4  =  0:  FLAGS  =  0 
40  PRINT  :  INPUT  "submarine  speed  in  knots";  V 

50  PRINT  :  INPUT  "distance  between  the  end  points  in  nautical  miles";  0 
■60  PRINT  ;  INPUT  ''track  length  in  nautical  miles";  TL 
70  IF  TL  <  0  THEN  60 

80  DN  «  SQR(TL  *  TL  *  D  *  0)  /  2:  REN  maximum  distance  from  the  x-axis 

90  SL  «  ON  /  7;  IF  0  /  2  >  SL  THEN  SL  »  (D  /  2)  /  INT(0  /  2  /  SL):  REN  length  of  a  cell  side 

100  TS  >  SL  /  V:  REN  time  step  in  hours 

110. NNS  •  30:  REN  imxi mum  number  of  steps 

120  NS  •  CINT(TL  /  SL):  REN  milter  of  time  steps  in  the  track 

130  IF  NS  >  NNS  THEN  GOTO  60 

140  SF  »  (64  /  35)  *  ON  /  1000:  REN  screen  factor 

150  AS .»  "cell  dimension  »  "  +  STRSfSL)  +  "  nm";  BS  «  "time  step  *"  STRS(TS)  *  "  hr" 

160  PRINT  -:  PRINT'  AS:  PRINT  :  PRINT  BS 

170  PRINT  :  INPUT  "delta  in  nautical  miles";  DEL;  REN  maximun  average  leg  length  in  nautical  miles 

180  IF  DEL  <  TL  •  TL  /  (TL  -  0)  /  50  THEN  170 

190  PRINT  ;  INPUT  "leg  length  distribution  index";  NI 

200  AS  »  "  PRINT  ;  INPUT  "standard  or  auxiliary  random  nunter  generator  (s/a)";  AS 

210  IF  AS  >  "A"  OR  AS  *  "a"  THEN  GOSUB  7720:  GOTO  280 

220  IF  AS  «  "S"  OR  AS  <  "s"  THEN  230  ELSE  200 

230. AS  a  "  ":  PRINT  :  INPUT  "supply  a  random  nunter  seed  (y/n)";  AS 

240  IF  AS  a  "N"  OR  AS  a  »n"  THEN  280 

250  IF  AS  a  "V'.'  OR- AS  a  Hyi.  thEN  FLAG4  a  1  ELSE  230 

260  PRINT  :  INPUT  "random  nunter  seed";  RNS 

270  IF  RNS  <  -327681  OR  RNS  >  32767  THEN  260  ELSE  RANOONIZE  RNS 

280  AS  a  H  H;  PRINT  :  INPUT  "view  only  single  tracks  (y/n)";  AS 

290  IF  AS  a  "YH  OR  AS  a  lyi  THEN  FLAG5  a  1;  NT  a  1;  GOTO  320 

300  IF  AS  a  "N"  OR  AS  a  Moh  THEN  310  ELSE  280 

310  PRINT  ;  INPUT  "nuiter  of  tracks";  NT 

320  AO  a  tL  /  2:  REN  semimajor  axis  of  the  initial  focusing  ellipse 

330  CO  a  0  /  2:  REN  distance  from  the  center  to  a  focus  of  the  initial  focusing  ellipse 

340  NNL  a  30 

350  LL  a  9:  LU  a  9:  NL  a  10:  NU  a  10:  NT  a  6 
360  REDIN  NP(-LL  TO  LU,  -NL  TO  NU,  HNS) 

370  REDIN  NT(-1  TO  1,  -1  TO  1,  -LL  TO  LU,  -NL  TO  HU,  NNS) 

380  REDIN  HTS(-1  TO  1,  -1  TO  1,  -LL  TO  LU,  -NL  TO  HU) 

390  FOR  J  a  1  TO  NT 

400  REDIN  X(HNL  ♦  1),  Y(MNL  ♦  1),  XS(MNL  +  1),  YS(MNL  ♦  1),  FR(MNL  ♦  1),  T(HNL  *  1),  TLEG(HHL  +  1) 

410  X(0)  a  -D  /  2:  Y(0)  a  0:  REN  coordinates  of  the  initial  point  for  all  tracks 

420  FR(0)  a  0:  TLEG(O)  a  o 

430  FLAG1  a  0 

440  A  a  AO:  C  a  CO:  S  a  0 

450  GOSUB  4770 

460  FOR  I  a  1  TO  MNL:  REN  leg  nunter 
470  B1  a  S0R((TL  +  0)  /  (TL  •  D)) 

480  82  a  SQR(1  -  D  *  D  /  TL  /  TL) 

490  GOSUB  7700 

500  ON  ERROR  GOTO  510:  GOTO  520 

510  RESUHE  490 

520  FLAG2  >  0 

530  IF  RAND  <a  .5  THEN  550 

540  RAND  a  RAND  -  .5:  FLAG2  a  1 

550  RC  a  2  *  ATN{B2  *  TAN(PI  *  RAND  -  ATN(BI))  +  D  /  TL):  REM  course  in  radians  relative  to  the  minor  axis  of 

the  focusing  ellipse 

560  ON  ERROR  GOTO  0 

570  IF  FLAG2  a  1  THEN  RC  a  pi  .  RC 

580  AL  a  (TL  -  0)  /  (TL  -  0  *  SIN(RC)):  :  REM  average  leg  length  for  the  relative  course  RC 
590  GOSUB  7620 
600  ON  ERROR  GOTO  0 

610  FR(I)  a  RC  +  S:  REH  course  in  radians  on  leg  I 

620  TLEGM  a  (A  *  A  -  C  *  C)  /  (A  -  C  *  SIN{RC)):  REH  maximum  leg  length  in  nautical  miles  for  the  relative 
course  RC 

630  IF  TLEG  >a  TLEGM  THEN  TLEG  =  TLEGM;  FLAG1  =  1:  REM  track  is  terminated  at  the  final  point 
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640-g6sUB' ^30 

65d‘Li>iEv,(X{ro,  - 

660jA  REM  r.eniainiog'trMk  length^-in  nautical  miles 

670  C  s'CI:  „REH:dist8tKe7tetw<(en  the  final  point  a^  the  end  point  of  the. previous  leg 
680’S  e  BR  -  PI'  /  2:  REM  direction  of  the  miiior  axis  of  the  focusing  ellipse 
:6904LEG(i')  »  TLEG;  REM  leg  length. of -leg  i'-' 

,700  IF  FLAG1  i  1  THEN  750 
-710.  NEXT  I 
:720.'NfR:«NTR +' 1 

730  .j;.  J  -  -1 

740' GOTO; 1280: ;REM<track  terminated 
TSO.'LINE  (0  /  2,..0)-.(X(I),  Y<I)) 

760'NLEG:»  I  1:  REM  ninber  of  legs 
<770  TlEG,«  2  *  C,:  TLECINLEG).  «  TLEG 

780  H  «  S  ♦  PI  /  2;"REM  bearing  in, radians  of  final  point  from  the  end  point  of  the  leg 

790  FR(NLEG)  «  H:'REM  course  in  radians  ;on  last  leg 

800  X(NIeG).«'0  /  2:  REM  x^cobrdinate  of  final  point  for  alt  tracks 

810  YCNLEG), »  0;:  REM  y-coordinate.df  final  point  for  all  tracks 

820  T(0)  *  0 

830  FOR  I  »  1 -TO  NLEG 

840  T,»  TLEGd)  /  V;  REM  time  in  hours  to  traverse  leg,  I 

850'T(I)  *  T(I  •  1)  ♦  T;  REM  time  in  hours  from  the  initial  point  to  the  end  point  of  leg  I 
860  NEXT  I 

870  X8(0)  *  -D  /  2;  YS(0)  *  0 
880  K  «,1;  TSK  »  TS 
890  FOR  I- »  1  TO  NtEG 

900  IF-  TSK  >  T(l)  AND  I  <  NLEG  THEN  960:  REM  the  next  time  step  point  is  past  the  Ith  leg 

910  TK.*  TSK  •  t(I  -  1):  REM  time  from  the  starting  point  of  the  Ith  leg  to  the.Kth  time  step  point 

920  XS(K)  »  X(I  -  1)  ♦  V  *  TK  *  S1N(FR(1))|  REM  x-coordinate  of  the  Kth  time  step  point 

930  YS(K)  *  Y{I  •  1)  ♦  V  *  TK  *  C0S{FR(I)):  REM  y-coordinate  of  the  Kth  time  step  point 

940. K-=  K  +  1:  TSK  *  TSK  t  TS:  REM  next  time  step 

950  IF  K  >  NiTHEN  980  ELSE  900 

960. NEXT  1 

970  ERASE  X,  Y,  FR;  T,  TLEG 
980  , FOR  K  »  1  TO  NS  -  1 
990’‘CIRCLE  <XS(K),  YS{K)),  8  •  SF 
1000  NEXT  K 

1010  IF  FLAGS  «  0  THEN  1050 

1020  AS  *  '•  INPUT  "Quit  (y/n)";  AS 

1030  IF  AS  «  "N"  OR  AS  »  "n"  THEN  CLS  :  GOTO  400 

1040  IF  AS  «  "Y"  OR  AS  »  "y"  THEN  END  ELSE  1020 

1050. FOR  1  *  0  TO  NS 

1060  FOR  L  *  -LL  TO  LU 

1070  IF  (L  •  1  /  2)  *  SL  <  XS(I)  AND  XS{1)  <=  <L  +  1  /  2)  *  SL  THEN  1090:  REM  the  x-coordinate  of  the  Ith  time 
Step  point  is  in  range  L 
1080  NEXT  L 

1090 -FOR  M  =  -ML  TO  MU 

1100  IF  (H  -  1  /  2)  *  SL  <  YS(I)  AND  YS(!)  <=  <M  ♦  1  /  2)  »  SL  THEM  1130:  REM  the  y-coordinate  of  the  Ith  time 
step  point  is  in  range  M 
1110  NEXT  M 
1120  GOTO  1270 

1130  MP(L,  M,  I)  *  MP(L,  M,  I)  +  1:  REM  random  tour  probability  map  element 

1140  IF  !  >  0  THEN  1170 

1150  LO  s  L:  MO  *  M 

1160  GOTO  1270 

1170  LI  »  L:  Ml  =  M 

1180  FOR  M  »  -1  TO  1 

1190  FOR  K  *  -1  TO  1 

1200  IF  NOT  (LI  =  LO  ♦  N  AND  Ml  =  MO  +  K)  THEN  1240 

1210  MT(N,  K,  LO,  MO,  I  -  1)  =  MT(N,  K,  LO,  NO,  I  -  1)  ♦  1:  REM  dynamic  diffusion  matrix  element 

1220  MTS(N,  K,  LO,  MO)  =  MTSfN,  K,  LO,  HQ)  +  1:  REM  static  diffusion  matrix  element 

1230  GOTO  1260 

1240  NEXT  K 

1250  NEXT  N 

1260  LO  »  LI:  HO  =  HI 

1270  NEXT  1 

1280  NEXT  J 

1290  AS  s  "  ":  PRINT  :  INPUT  "display  a  random  tour  probability  map  (y/n)";  AS 
1300  IF  AS  =  "N"  OR  AS  =  "n"  THEN  1500 
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1310  fF  M  AS  =  "y«  THEM  1320  ELSE  1290 

1320  SCREEN  0:  CLS 

.13M  PRINT  "MS  a  NS  . 

-13iO-PRiNT  :  INPUT  “lehter  the  number  of,  time  steps  from, 0  to  NS!';  I 
i350  IF  I'.?<  6  OR  T  >  NS  THEM  ISW 
13M'FOR,M  a  10  TO  -ilO  STEP  -1 

1370  X$-»  STRS(MP<-9,  M,  I));  X$  aLTRIM$(RTRIMS(Xi)):  PRINT  X$;  TAB(5); 

13M  FOR^N  a  ^8  TO  8 

1390  X$:«  STR$(MP(N,  M,  I»;  X$  a  LTRIM$<RTRIM$(X$));  PRINT  X$;  TAB((N  +  9)  •  4  +  5); 

1400  NEXT  N 

1410  X$  a  sfR$(MP(9,.M,.  D):  XS  aafRIM$(RTRIM$(X$)):  PRINT  XS 
142q.NEXT  M 

■1430  AS  a  "  PRINT  ;  INPUT  "print  the  random  tour  probability  map  (y/n)";  AS 

1440:IF  AS  i  «N''  OR  AS  a  Mnii  THEM  1290 

1450  IF  AS -a  «Y"  OR.  AS  a  , My.;  then  1460  ELSE  1430 

1460:WIDTH  LPRINT  130 

1470.GOSUB  4880 

1480  WIDTH  LPRINT  80 

1490  GOTO 

1500  AS  a  PRINT  ;  INPUT-  "print  a  random  tour  probability  map  (y/n)";  AS 

1510: IF  AS  a  mnh  or  as  a  MnM  JHEN  1610 

1520  IF  AS  a  MY"  OR  AS  «  wy"  THEN  1530  ELSE  1500 

1530  SCREEN  0:  CLS 

1540  PRINT  "NS  a  ";  NS 

1550  P'RINT  :  INPUT  "enter  the  nuiber  of  time  steps  from  0  to  NS  ";  I 

1560  IF  I  <  0  OR  I  >  NS  THEN  1550 

1570  WIDTH, LPRINT  130 

1580  COSUB  4880 

1590  WIDTH  LPRINT  80 

1600  GOTO  1290 

1610; AS  a  "  "s  PRINT  :  INPUT  "print  all  the  random  tour  probability  maps  (y/n)";  AS 

1620  IF  AS  a  "N"  OR  AS  a  "n"  THEN  1680 

1630  IF  AS  a  "Y"  OR  AS  a  "y"  THEN  1640  ELSE  1610 

1640-UIDTH  LPRINT  130 

1650  FOR  I  a  0  TO  NS 

1660  GOSUB  4880 

1670  NEXT  I 

1680  AS  a  "  ":  PRINT  ;  INPUT  "print  all  the  random  tour  probability  maps  to  a  file  (y/n)";  AS 

1690  IF  AS  a  "N"  OR  AS  a  "n"  THEN  1800 

1700  IF  AS  a  "Y"  OR  AS  a  "y"  then  1710  ELSE  1680 

17.10  PRINT  :  INPUT  "input  the  data  file  name";  FS 

1720  ON  ERROR  GOTO  1730:  GOTO  1740 

1730  RESUME  1710 

1740  OPEN  "0",  #1,  FS 

1750  FOR  I  a  0  TO  NS 

1760  GOSUB  7050 

1770  NEXT  I 

1780  CLOSE  #1 

1790  ON  ERROR  GOTO  0 

1800  ERASE  HP 

1810  AS  a  H  ":  PRINT  :  INPUT  "display  a  dynamic  transition  matrix  (y/n)";  AS 

1820  IF  AS  a  "N"  OR  AS  a  "n"  THEN  1920 

1830  IF  AS  a  "Y"  OR  AS  a  »y"  THEN  1840  ELSE  1800 

1840  SCREEN  0:  CLS  :  GOSUB  6480 

1850  PRINT  :  PRINT 

1860  FOR  K  a  1  TO  -1  STEP  -1 

1870  PRINT  ;  PRINT  HT(-1,  X,  L,  M,  I  •  1);  TAB(15);  MT(0,  X,  L,  H,  1  •  1);  TAB(29);  MT(1,  X,  L,  M,  1-1) 
1880  NEXT  X 

1890  PRINT  :  PRINT  :  PRINT  :  PRINT  "nunber  of  tracks  completed  =  ";  NT 
1900  PRINT  :  PRINT  "nunber  of  tracks  terminated  =  ";  NTR 
1910  PRINT  :  GOTO  1810 

1920  AS  a  "  PRINT  ;  INPUT  "print  a  dynamic  transition  matrix  (y/n)";  AS 

1930  If  AS  a  "N"  OR  AS  a  "n»  THEN  2100 

1940  IF  AS  a  "Y"  OR  AS  a  "y"  then  1950  ELSE  1920 

1950  SCREEN  0:  CLS  :  GOSUB  6480 

1960  AS  a  NS  ♦  "  dynamic  transition  matrix";  LPRINT  AS 
1970  LPRINT  ;  LPRINT  "transition  number  =  ";  I 
1980  LPRINT  ;  LPRINT  "x  =  ";  L 
1990  LPRINT  "y  a  ";  m 
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,-2bbO/F6R,jC  s;,i  fo  -i  STEP  -1  . 

;'2010;iPRlNT  LPRINT  ;  tAB(i5>r  Mf(-1,  K,.  L,  M,  I  -  1);  tAB(29);  MT(0,  K,  L,  M,  I  -  T);  TAB(43);  MT(1,  K,  L, 
M,„i:  r  i) 

'2020  NEXT 

203b  LPRINT:  LPRINT 
,  2040  GOSUB' 6130 

2050'JIF  FLAC4  =vO:THEN  N-=  34  ELSE  N  =  3*3 
:206b  FOR  K'S  I  TO  N  - 

.2070  LPRINT 
2080', NEXT  if  , 

■20W  bdfb  1920 

2100  A%  PRINT  :  INWT  Pprint  nonzero. dynamic  transi.tion  matrices  (y/n)*';  AS 

2110  IF  AS  «  «N'.'  OR  AS  *  "n"  THEN  2390 

2120  IF  AS  »  "Y''-OR  AS  •  "y”  THEN  2130  ELSE  2100 

2130  SCREEN  0:  CLS 

2140  PRINT  "NS  NS 

2150.  PRINT  :  INPUT  "enter  the  transition  nliiter  from  1  to  NS";  I 
216blF  I  <  1  OR  I  >  NS  THEN  2150 

2170  AS  »  NS  t  ".nonzero  dynamic  transition  matrices":  LPRINT  AS 
2180-]iPRTNT  :  LPRINT  "transition  number  =  ";  1:  LPRINT 
2190  FOR  L  »  -(LL  -  1)  TO  (LU,  -  1) 

2200  FOR  H  «  -{ML  -  1)  TO  (MU  -  1) 

2210, MAT  =  0 

2220  FOR  N  =  -1  TO  1 

2230  FOR  K  »  -1  TO  1 

2240  HAT  *  MT(N,  K,  L,  M,  I  -  1)  +  HAT 

2250  NEXT  K 

2260 'NEXT  N 

.2270  IF  MAT  *  0  THEN  2340 
■2280  LPRINT  :  LPRINT  "x  =  »;  L 
2290  LPRINT  "y  *  M 
2300  FOR  K  *  1  TO  -1  STEP  -1 

2310  LPRINT  ;  LPRINT  ;  TAB(15);  MT(-1,  K,  L,  M,  I  -  1);  TAB(29);  MT{0,  K,  L,  H,  I  -  1);  TAB(43);  MT(1,  K,  L, 
M,  I  •  1) 

2320  NEXT  K 

2330  LPRINT  :  LPRINT 

2340  NEXT  M 

2350  NEXT  L 

2360  LPRINT  :  LPRINT 

2370  GOSUB  6130 

2380  GOTO  2100 

2390  AS  »  "  PRINT  :  INPUT  "print  a  dynamic  transition  matrix  to  a  file  (y/n)";  AS 

2400  IF  AS  «  "N"  OR  AS  =  "n»  THEN  2630 

2410  IF  AS  =  "Y»  OR  AS  =  "y"  THEM  2420  ELSE  2390 

2420  SCREEN  0;  CLS  ;  INPUT  "input  the  data  file  name";  FS 

2430  ON  ERROR  GOTO  2440:  GOTO  2450 

2440  RESUME  2420 

2450  OPEN  "0",  #1,  FS 

2460  GOSUB  6480 

2470  AS  =  NS  "  dynamic  transition  matrix":  PRINT  #1,  AS 
2480  PRINT  #1,  :  PRINT  #1,  "transition  number  =  ";  I 
2490  PRINT  #1,  :  PRINT  #1,  "x  *  ";  L 
2500  PRINT  #1,  "y  =  ";  M 
2510  FOR  K  «  1  TO  -1  STEP  -1 

2520  PRINT  #1,  :  PRINT  #1,  ;  TAB(15);  HT(-1,  K,  L,  M,  I  -  1);  TAB(29);  HT(0,  K,  L,  M,  I  -  1);  TAB(43);  MT(1, 
K,  L,  M,  I  -  1) 

2530  NEXT  K 

25 vO  PRINT  #1,  :  PRINT  #1, 

2550  GOSUB  6560 

2560  IF  FLAG4  =  0  THEN  N  =  27  ELSE  N  =  26 
2570  FOR  K  *  1  TO  H 
2580  PRINT  «1, 

2590  NEXT  K 
2600  CLOSE  #1 
2610  ON  ERROR  GOTO  0 
2620  GOTO  2390 

2630  AS  =  "  PRINT  :  INPUT  "generate  dynamic  diffusion  tracks  (y/n)";  AS 

2640  IF  AS  *  "Y"  OR  AS  =  "y"  THEN  2660 

2650  IF  AS  =  "N"  OR  AS  =  "n"  THEN  2770  ELSE  2630 
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■26W):%CREEN;;6:'{CLS; 

2i70>|RINj  J'iis  l|S 

^6&l)1>RINf  INWiX '('enter  the- nuriijer.of  diffusions  from  1  to  NS";  ND 
i26?q:dFVNbv<  i;dR;N0;>^;N^ 

.2700; PRINT  INPUT  "enter  the  hunter  of  trecks";  NOT 

;2710CLS  ,  .  .  - 

2720  ObSUB  4770 

.2730  FOR.  J,  «  1  to  NOT 

■2740  GOSUB  5050 

2750  riEXT  J  . 

2760  OOTO  2«0 

2770  A*  »  "  •<;  PRINT  i  INPUT,  generate dynamic  diffusion  probability  map  (y/n)";  AS 

2780  tF  . AS  ■.47"  OR  AS  ■  »y"  THEN  2800 

2790.  IF  AS  »  «N"  OR  AS  «  "n"  THEN  3050  ELSE  2770 

280Ci;SCREEN  0:  -CLS 

28)0  PRINT  "NS  »  NS 

2820  PRINT  :  INWT''!enter  the  huiter  oif  diffusions  from  0  to  NS  ND 

2830  IF  NO  <  0  OR  NO  >  NS  THEN  2820 

2840  REDIM  NPO(-LL  TO/LU,  -ML  TO  HU)  AS  INTEGER 

2850  GOSUB  5330 

2860  AS'«  "  PRINT  ;  INPUT  "display  the  dynamic  diffusion  probability  map  (y/n)";  AS 

2d70'ir  AS  »  "N"  OR  AS-»  «h"  THEN  2970 

2880  IF  AS  »  «Y"  OR  AS  «  "y"  THEN  2890  ELSE  2860 

2890  SCREEH^O:  CLS 

2900  FOR  H  «  10  TO  -10  STEP  -1 

'2910;XS  »  STRS(HP0<'9,  M));  XS  «  LTRIMS(RTRIMS(XS)):  PRINT  XS;  TAB(5); 

2920  FOR-N  »  '8  TO  8 

2930  XS4  STRS(MP0(N,  M)):  XS  «  LTRIHS(RTRiMS(XS}):  PRINT  XS;  TAB((N  ^  9)  •  4  ^  5); 

2940' NEXT  N 

2950.  XS’«  STRS(NP0(9,  M)):  XS  «  LTRIMS(RTRIKS(XS)):  PRINT  XS 
2960.NEXT  H 

2970.' ^S  «•"  PRINT  :  INPUT  "print  the  dynamic  diffusion  probability  map  (y/n)";  AS 

,2980  IF  AS  »  "N"  OR  AS  »  "n"  THEN  2770 

2990;'1F  AS'«  "Y"  OR  AS  «  "y"  THEN  3000  ELSE  2970 

3000  WIDTH' LPRINT  130 

301O'COSUB  5680 

3020  WIDTH  LPRINT  80 

3030  ERASE  MPO 

3040. GOTO  2770 

3050  AS  »  "  "'.  PRINT  :  INPUT  "print  all  the  dynamic  diffusion  probability  maps  (y/n)";  AS 

3060  IF  AS  «  "N"  OR  AS  »  «n"  THEN  3160 

3070  IF  AS'»  "Y«  OR  AS  »  "/•  THEN  3080  ELSE  3050 

3080  WIDTH  LPRINT  130 

3090  FOR  ND  «  0  TO  NS 

3100  REOIH  MP0(-LL  TO  LU,  -ML  TO  HU)  AS  INTEGER 

3110  GOSUB  5330 

3120.C0SU8  5680 

3130  NEXT  NO 

3140  ERASE  MPO 

3150  WIDTH  LPRINT  80 

3160  AS  »  "  PRINT  ;  INPUT  "print  all  the  dynamic  oiffusion  probability  maps  to  a  file  (y/n)";  AS 

3170  IF  AS  *  "N»  OR  AS  »  "n"  THEN  3310 

3180  IF  AS  »  »Y"  OR  AS  *  «y"  THEN  3190  ELSE  3160 

3190  PRINT  :  INPUT  "input  the  data  file  name";  FS 

3200  ON  ERROR  GOTO  3220 

3210  GOTO  3230 

3220  RFSI.'ME  3190 

3230  ^PEM  "0",  #1,  fS 

324C  FCR  NO  3  0  TO  NS 

3250  REOIM  MPDCLL  TO  LU.  -ML  TO  MU)  AS  INTEGER 

3260  uOSUB  5330 

3270  GOSUB  6700 

3280  NEXT  NO 

3290  CLOSE  «1 

3300  ON  ERROR  GOTO  0 

3310  ERASE  MPD 

3320  AS  =  "  ":  PRINT  :  INPUT  "display  a  static  transition  matrix  (y/n)";  AS 

333P  IF  AS  =  "N"  OR  AS  =  "n"  THEN  3430 

3340  If  AS  =  "Y"  OR  AS  =  "y"  THEN  3350  ELSE  3320 
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:3356-SCREEN:.0:  CLS  ; 'GOSUB'6510: 

3360 -  PR  I  NT  -;  .PR  I  NT. 

3370  FOR  1C  *  1  TO  -1  STEP  -1 

;3380..PR1NT  ;  PRINT  MtS(,-1,  K,  L,  M>;  TAB(15);  MTS(0,  K,  L,  M);  TAB(29);  MTS(1,  K,  L,  H) 

3390.  NEXt  N 

3400  PRINT  :  PRINT,  ;  PRINT  :  PRINT  "nurtber  of  tracks  completed  =  NT 
34l0  PRINT  ;  PRINT  "number  of  tracks  terminated  =  NTR 

3420  PRINT-  :  GOTO-; 3320 

3430  A$  «.  ••  PRINT  ;  INPUT  "print  a. static  transition  matrix  (y/n)";  AS 

3440' IF  AS  »  «N"  OR  AS  »'«n"  THEN  3600 

3450  IF  AS  »  "Y"  OR  AS  «  "y"  THEN  3460  ELSE  3430 

3460  SCREEN  0:  CLS  GOSUB  6510 

3470  AS  «  NS  +  "  static  transition  matrix";  LPRINT  AS 

3480  LPRINT  ;  LPRINT  "x  »  L 

3490  LPRINT  «y  »  H 

3500  FOR  K  ■  1  TO  -1  STEP  -1 

3510  LPRINT  ;  LPRINT  ;  TAB(15);  NTS(-1,  K,  L,  M);  TAB(29);  MTS{0,  K,  L,  M>;  TAB(43);  MTS(1,  K,  L,  M) 

3520  NEXT  K 

3530  LPRINT  ;  LPRINT 

3540  GOSUI  6130 

3550  IF  FLAG4  ■  0  THEN  N  >  36  ELSE  N  ■  35 

3560  FOR  1C  »  1  TO  N 

3570  LPRINT 

3580  NEXT  K 

3590  GOTO  3430 

3600  AS  »  "  «;  PRINT  ;  INPUT  "print  all  noniero  static  transition  matrices  (y/n)";  AS 

3610  IF,  AS  ■  "N"  OR  AS  »  "n»  THEN  3840 

3620  IF  AS  «  "Y"  OR  AS  «  "y"  THEN  3630  ELSE  3600 

3630  SCREEN  0;  CLS 

3640  AS  «  NS  ♦  "  noniero  static  transition  matrices";  LPRINT  AS 
3650  FOR  L  »  -(LL  •  1)  TO  (LU  •  1) 

3660  FOR  M  •  -(HL  -  1)  TO  (MU  •  1) 

3670  MAT  •  0 

3680  FOR  N  ■  -1  TO  1 

3690  FOR  K  «  *1  TO  1 

3700  MAT  >  MTS(N,  K,  L,  M)  t  MAT 

3710  NEXT  K 

3720  NEXT  N 

3730  IF  MAT  «  0  THEN  3800 
3740  LPRINT  ;  LPRINT  "x  »  L 
3750  LPRINT  "y  «  »;  M 
3760  FOR  K  »  1  TO  -1  STEP  *1 

3770  LPRINT  ;  LPRINT  ;  TAB(15);  MTS(-1,  K,  L,  M);  TAB(29);  MTS(0,  K,  L,  M);  TAB(43);  MTS(1,  X,  L,  M) 

3780  NEXT  K 

3790  LPRINT  ;  LPRINT 

3800  NEXT  M 

3810  NEXT  L 

3820  LPRINT  ;  LPRINT 

3830  GOSUB  6130 

3840  AS  »  "  ";  PRINT  ;  INPUT  "print  a  static  transition  matrix  to  a  file  (y/n)";  AS 

3850  IF  AS  »  «N"  OR  AS  *  "n»  THEN  4080 

3860  IF  AS  »  "Y"  OR  AS  =  "y"  THEN  3870  ELSE  3840 

3870  SCREEN  0;  CLS  ;  INPUT  "input  the  data  file  name";  FS 

3880  ON  ERROR  GOTO  3890;  GOTO  3900 

3890  RESUME  3870 

3900  OPEN  "0",  #1,  FS 

3910  GOSUB  6510 

3920  AS  3  NS  -f  "  static  transition  matrix" 

3930  PRINT  «1,  AS 

3940  PRINT  #1,  ;  PRINT  #1,  "x  =  ";  L 
3950  PRINT  #1,  "y  »  ";  M 
3960  FOR  X  »  1  TO  -1  STEP  -1 

3970  PRINT  #1,  ;  PRINT  #1,  ;  TAB(15);  MTS(-1.  X,  L,  M);  TAB(29);  MTS(0.  X,  L,  M);  TAB(43);  MTS(1,  X,  L,  M) 
3980  NEXT  K 

3990  PRINT  #1,  ;  PRINT  #1, 

4000  GOSUB  6560 

4010  IF  FLAG4  =  0  THEN  N  =  30  ELSE  N  =  29 
4020  FOR  X  =  1  TO  N 
4030  PRINT  #1, 
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4050  CLOSE  #1 
40M;ON  , ERROR  GOTO' 5 
4070.  GOTO  ;3856 

4080  AS  "  T*  PRINt  :  INPUT  "generate,  static  diffusion  tracks  (y/n)";  AS 

4090  IF.  AS  »  "Y-"  OR  AS  ?  "y«'THEN'4110 

4100  IF  AS.s  "N«=OR  AS  i  "n'.'  THEN  4230  ELSE  4080 

4110  SCREEN  O':  CLS 

4120  PRINT  "NS  «  NS 

4130  PRINT  :  INPUT  "enter  the  nuiter  of  diffusions  from  1  to  NS";  NO 

4140  IF  NO  <  1  OR  NO  >  NS  THEN  4130 

4.150  I^RINT  :  INPUT  "enter  the  ntmber  of  tracks";  NOT 

4160  CLS 

4170  ’ 4770 

4180  FLAG3  >  1 

4190  FOR  J  »  1  TO  NOT 

4200  GO^B  5050 

4210  NEXT  J 

4220  GOTO  4080 

<4230  AS  *  "  PRINT  :  INPUT  "generate  a  static  diffusion  probability  map  (y/n)";  AS 

4240  IF  AS  »  "Yi*  OR  AS  »  "y"  THEN  4260 

4250  IF  AS  «  "NH  OR  AS  «  «n"  THEN  4520  ELSE  4230 

4260  SCREEN  0:  CLS 

4270;PRINT  "NS  «  ";  NS 

4280  PRINT' :  INPUT  "enter  the  number  of  diffusions  from  0  to  NS  ";  NO 
4290  IF  NO  <  0  OR  NO  >  NS  THEN  4280 
4300‘RE0IH  P0T<'LL  TO  LU,  -ML  TO  HU) 

4310  GOSUB  7220 
4320  GOSUB  7550 

4330  AS  «  "  ":  PRINT  :  INPUT  "display  the  static  diffusion  probability  map  (y/n)";  AS 

4340  IF  AS  «  "N"  OR  AS  >  "n"  THEN  4440 

4350  IF  AS  >  "Y"  OR  AS  «  "y"  THEN  4360  ELSE  4330 

4360  SCREEN..0:  CLS 

4370  FOR  M  •  10  TO  -10  STEP  -1 

4380  XS  «  $TRS(CINT(P0T(-9,  M))):  XS  «  LTRIMS(RTRIHS(XS)):  PRINT  XS;  TAB(5); 

4390  FOR  N  ■  *8  TO  8 

4400  XS  •  STRS(CINT(POT(N,  M»):  XS  «  LTRIMS(RTRIMS(XS)):  PRINT  XS;  TAB((H  ♦  9)  *  4  +  5); 

4410  NEXT  N 

4420  XS  «  STRS(CINT(P0T(9,  H))}:  XS  «  LTRIHS(RTRIHS(XS}):  PRINT  XS 
4430  NEXT  H 

4440  AS  »  "  ":  PRINT  ;  INPUT  "print  the  static  diffusion  probability  map  (y/n)";  AS 

4450  IF  AS  ■  »N"  OR  AS  »  "n»  THEN  4230 

4460  IF  AS  *  "Y"  OR  AS  «  "y"  THEN  4470  ELSE  4440 

4470  WIDTH  LPRINT  130 

4480  GOSUB  5850 

4490  WIDTH  LPRINT  80 

4500  ERASE  PDT 

4510  GOTO  4230 

4520  AS  *  "  ":  PRINT  ;  INPUT  "print  all  the  static  diffusion  probrbility  maps  (y/n)";  AS 

4530  IF  AS  »  "N"  OR  AS  »  "n"  THEN  4630 

4540  IF  AS  »  "Y"  OR  AS  *  "y"  THEN  4550  ELSE  4520 

4550  WIDTH  LPRINT  130 

4560  FOR  ND  >  0  TO  NS 

4570  REDIH  P0T(*LL  TO  LU,  -ML  TO  MU) 

4580  GOSUB  7220 
4590  GOSUB  5850 
4600  NEXT  NO 
4610  ERASE  PDT 
4620  WIDTH  LPRINT  80 

4630  AS  *  "  ":  PRINT  ;  INPUT  "print  all  the  static  diffusion  probability  maps  to  a  file  (y/n)" 

4640  IF  AS  s  "N"  OR  AS  s  "n"  THEN  4760 

4650  IF  AS  s  "Y"  OR  AS  =  "y"  THEN  4660  ELSE  4630 

4660  PRINT  :  INPUT  "input  the  data  file  name";  FS 

4670  ON  ERROR  GOTO  4680:  GOTO  4690 

4680  RESUME  4660 

4690  OPEN  "0",  #1,  FS 

4700  FOR  ND  =  0  TO  NS 

4710  REDIM  POT(-LL  TO  LU,  -ML  TO  MU) 

4720  GOSUB  7220 


-4730: COSUB :6870 
4740- NEXT  NO  ' 

4750-.ERASE  POT 
'4760.'END 

4770  XM  «  iO  •  SLj  YH  =  8  *  SL 
4780  screen: 12 

4790‘.WiN0OW  (-XM,  -YM)-<XM,  YM) 

4800  CIRCLE  (-0  /  2,  0),  10  *  SF;  CIRCLE  <0/2,  0),  10  *  SF 
4810  FOR  M  *  0  TO  19 

4'82d;LINE  C-XM  +  (M  +  1  /  2)  *  SL,  YM  -  SL  /  2)-(-XM  ♦  (M  ♦  1  /  2)  *  SL,  -YM  +  SL  /  2) 

4830  NEXT  M 

4840  FOR  M, «  0  TO  15 

4850  LINE  (-XM  ♦  SL  /  2,  YM  •  <M  +  1  /  2)  •  SL)-(XM  -  SL  /  2,  YM  •  (M  +  1  /  2)  *  SL) 

4860  NEXT  M 
4870  RETURN 

4880’AS  »  NS  ♦  random  tour  probability  map":  LPRINT  AS 
4890  LPRINT  ;  LPRINT  "nunber  of  tima  steps  *  I:  LPRINT 
4900  FOR  M  «  10  TO  -10  STEP  -1 

-4910  XS  »  STRSCMP(-9,  M,  I));  XS  s  LTRIMS(RTRIMS(XS)):  LPRINT  XS;  TAB(MT); 

4920  FOR  N  '*  -8  TO  8 

4930  XS  »  STRS(MP(N,  M,  I»;  XS  «  LTRIMS<RTRIMS(XS)>;  LPRINT  XS;  TAB{(N  ♦  9)  *  (MT  -  1)  +  MT) 
4940  NEXT  N 

4950  XS  S  STRS(NP(9,  M,  1)}:  LPRINT  XS 
4960  NEXT'  M 
4970  GOSUB  6130 
4900  COSUB  6270 

4990APRINT  "cell  entry  sum  »  SM 

5000  IF  FLAG4  *  0  THEN  N  *  22  ELSE  N  »  21 

5010- FOR  K  s  1  TO  N 

5020  LPRINT 

5030  NEXT  K 

5040.  RETURN 

5050  LO."  *0  /  2  /  SL;  HO  »  0;  REM  initial  cell  indices 
5060  FOR  I  »  1  TO  ND 

5070.RE0IM  MTA(-1  TO  1,  -1  TO  1)  AS  INTEGER 

5080  MPSUM  «  0 

5090 -FOR  N  »  -1  TO  1 

5100  FOR  X  *  -1  TO  1 

5110  IF  FLAG3  »  0  THEN  5130 

5120  MPSUM  s  MTS(N,  K,  LO,  MO)  +  MPSUM:  GOTO  5140 

5130  MPSUM  »  MT(N,  X,  LO,  MO,  I  •  1)  ♦  MPSUM 

5140  MTACN,  X)  s  MPSUM 

5150  NEXT  K 

5160  NEXT  N 

5170  COSUB  7700 

5180  RNU  >  RAND  *  MPSUM 

5190  FOR  N  *  -1  TO  1 

5200  FOR  X  *  -1  TO  1 

5210  IF  RNU  <  MTACN,  X)  THEN  5260;  REM  transition  from  (L0+N,M0+X)  to  (L0,M0) 

5220  NEXT  X 
5230  NEXT  N 
5240  ERASE  HTA 
5250  GOTO  5320 

5260  LI  s  LO  >  N:  Ml  s  MO  +  K 

5270  LINE  (LI  *  SL,  Ml  •  SL)-<L0  *  SL.  MO  *  SL) 

5280  CIRCLE  (LI  *  SL,  HI  *  SL),  8  *  SF 

5290  LO  s  LI;  MO  =  Ml 

5300  ERASE  MTA 

5310  NEXT  I 

5320  RETURN 

5330  REOIM  MSUM(-LL  TO  LU,  -ML  TO  MU,  MNS) 

5340  FOR  I  »  0  TO  NS  -  1 
5350  FOR  L  =  -LI  TO  LU 
5360  FOR  M  =  -ML  TO  MU 
5370  FOR  N  =  -1  TO  1 
5380  FOR  X  =  -1  TO  1 

5390  MSUHCL,  M,  I)  =  HT(N,  X,  L,  M,  I)  +  MSUH(L,  H,  I) 

5400  NEXT  X 
5410  NEXT  N 
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5420  NEXT  M 
5430;  NEXT  b; 

5440  NEXT  I 

5450  MroCrD  /  2  /  SL,  6)  *  NT:  REM  initiali^eelT  probability  factor 

5460  IF  N0  *  d  THEN  5660 

5470  FOR  I  X  1  TO  NO 

5480/REblM  Bl-LL  TO'  LU,  -ML  TO  MU), 

5490- FOR  Lx  -<LL  -  1)  TO  (LU  -  1) 

5500  FOR  M'x  .(ML  -  1)  TO  (MU  -  1) 

5510jFOR  N  X  -1  TO  1 
5520  FOR  K  X  -1  TO  1 

5530  IF  MSUM(L  •  M,  M  -  Kj  I  -  1)  *  6  THEN  5550 

5540  8(L,  M)  x.'MT(N,  K,  I  -  N,  M  ■  K,  !  -  1)  /  MSUM(L  -  N,  M  -  K,  I  -  1)  *  MP0(L  -  N,  M  -  K)  +  B(L,  M) 

5550  NEXT  K 

5560 - NEXT  N 

5570  NEXT  M 

5580  NEXT  L 

5590  F^  L  X  -LL  TO'LU 

5600  FOR  M  X  -ML  TO- MU 

.5610.^D(L,  M)  X  B(L,.M}:  REM  dynamic  diffusion  probability  map  element 

'5620;  NEXT  M 

5630  NEXT  L 

5640  ERASE  B 

5650  NEXT'  I 

5660 .ERASE  MSUM 

5670'.RETURN 

.5680  AS  x.NS  +  !'-dynamic  diffusion  probability  map";  LPRINT  AS 
5690  L'PRINT  :  LPRINT  'Inumber  of  diffusions  »  «;.nD:  LPRINT 
5700  FOR  M  X  10  TO  -10  STEP  -1 

5710  XS  X  STRS(MPb(-9,  M)):  XS  »  LTRIMS<RTRIMS(XS)):  LPRINT  XS;  TAB(KT); 

5720  FOR  N  X  -8  TO  8 

5730  XS  X  STRS(MP0(N,  M));  XS  x  LTRIMS(RTRIMS(XS)):  LPRINT  XS;  TAB((N  +  9)  *  (MT  •  1)  +  MT); 

5740  NEXT  N 

5750,  XS  X  STRS(MP0(9,  M)):  XS  x  LTRIMS(RTRIMS(XS)):  LPRINT  XS 
5760  NEXT  M 
5770  COSUB  6130 
5780  GOSUB  6340 

5790  LPRINT  "cell  entry  sum  «  SM 

5800  IF  FLAG4  x  Q  THEN  N  *  22  ELSE  N  x  21 

5810  FOR  K  X  1  TO  N 

5820  LPRINT 

5830  NEXT  K 

5840  RETURN 

5850  AS  X  NS  "  Static  diffusion  probability  map":  LPRINT  AS 
5860  LPRINT  :  LPRINT  "number  of  diffusions  «  ";  ND:  LPRINT 
5870  FOR  M  X  10  TO  -10  STEP  -1 

5880  XS  X  STRS(CINT(PDT(-9,  M))):  XS  x  LTRIMS(RTRIMS(XS)):  LPRINT  XS;  TAB(MT); 

5890  FOR  N  X  -8  TO  8 

5900  XS  X  STRS(CINT(POT(N,  M))):  XS  x  LTRIMS(RTRIMS(XS)):  LPRINT  XS;  TAB((N  +  9)  *  (MT  -  1)  +  MT); 

5910  NEXT  N 

5920  XS  X  STRS(CINT(P0T(9,  M))):  XS  x  LTRIMS(RTRINS(XS)):  LPRINT  XS 
5930  NEXT  M 
5940  GOSUB  6130 
5950  GOSUB  6410 

5960  LPRINT  "cell  entry  sum  x  sM 

5970  LPRINT  "cell  entries  are  rounded  integer" 

5980  IF  FLAG4  x  0  THEN  N  x  21  ELSE  N  x  20 

5990  FOR  K  X  1  TO  N 

6000  LPRINT 

6010  NEXT  K 

6020  RETURN 

6030  X(l)  X  x(I  -  1)  +  TLEG  *  SIN(FR(I)):  REM  x-coordinate  of  the  end  of  leg  I 

6040  Yd)  X  Y(1  -  1)  +  TLEG  *  COS(FR(I)):  REM  y-coordinate  of  the  end  of  leg  I 

6050  X  X  0  /  2  -  X(I):  Y  =  0  -  Y(I):  REM  determine  the  bearing  and  range  of  the  final  point  from  the  end  of  leg 
I 

6060  RF  X  SQR(X  *  X  +  Y  *  Y):  REM  range  of  the  final  point  from  the  end  of  leg  I 

6070  Cl  X  RF  /  2:  REM  distance  of  the  center  of  the  trial  focusing  ellipse  from  the  final  point 

6080  AT  X  A  -  TLEG  /  2;  REM  one-half  the  remaining  track  length 

6090  If  A8S(X  /  RF)  X  1  THEN  83  x  pj  /  2  *  SGN(X)  ELSE  88  =  ATN(X  /  RF  /  SQR(1  -  X  *  X  /  RF  /  RF)) 
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rilOO  IF-  ABS(Y  /  RF)/?  I’  THEN  BR  =  PI  /  ,2  *  (1  -  SGM{Y))  ELSE  BR  =  PI  /  2  -  ATNIY  /  RF  /  SOR(1  -  Y  *  Y  /  RF  / 

6H0,';IF  BB  <  0  THEN  BR  =  2  *  PI  -^BR;  REM  bearing  in  r^ians  of  the  final  point  from  the  end  point  of  leg  I, 
6j20'RETURN 

6130.‘LPRiNT', :  tPRINT.  "submarine  speed:  in  knots  =  'i;  V 

6140; tPRINT'  "distance- tetween: end  points  in  nautical  miles  =  ";  o; 

6150  L^RINT  "length  of  a  cell  side  in  nautical  miles  =  ";  SL 

.6160.  LPRINT  "time  step  in,  hours  ='!;  TS 

6170  LPRiNT' "maxiimiii  midbef  of  time-steps  =  ";  NS 

'6180  LPRINT  "track  length  in  nautical  miles  »  ";  TL 

6i90''LPRINT  "delta  in  nautical /miles  >  DEL 

;620b  (.PRINT  ".l«^j  length  distHibution  index  =  «;  NI 

.6210  LPRINT  Hrandom .nirtber  generator  «  «;  R$ 

6220  IF/FLAW  ■  O'THEN  6240 
6230 -LPRINT  Pr,ahdoffl  number  seed  «  ";  RNS 
6240  LPRINT  Hrxurber  of  tracks  completed  *  NT 
6250  LPRINT  "huiber  of  tracks  terminated  =  ";  NTR 
6260  RETURN 
6270  SM  »  0 
6280  FOR  L  ■  -LL  TO,  LU 
6290  FOR  M  «  -ML  TO  MU 
6300 <SM  «  MP<L,  M,  I)  ♦  SM 
6310  NEXT  M 
6320;  NEXT  L 
6330  RETURN 
•6340  SM  *0 
6350  FOR  L  *  -LL  TO  LU 
6360- FOR  M  *  -ML  TO  MU 
6370 'SM«  MP0(L,  M)  +  SM 
6380  NEXT  M 
6390  NEXT  L 
6400  RETURN 
6410  SM  »  0 
6420  FOR  L  »  -LL  TO  LU 
6430  FOR  M  »  -ML  TO  MU 
6440  SM  «  P0T(L,  M)  +  SM 
6450  NEXT  M 
6460  NEXT  L 
M70  RETURN 
6480  PRINT  "NS  *  NS 

6490  PRINT  ;  INPUT  "enter  the  transition  number  from  1  to  NS";  I 

6500  IF  I  <  1  OR  1  >  NS  THEN  6490 

6510  PRINT  5  INPUT  "x-coordinate";  L 

6520  IF  L  <  -ILL  -  1)  OR  L  >  (LU  •  1)  THEN  6510 

6530  PRINT  !  INPUT  "/-coordinate";  M 

6540  IF  M  <  -(ML  -  1)  OR  M  >  (MU  ■  1}  THEN  6530 

6550  RETURN 

6560  PRINT  #1,  ;  PRINT  #1,  "submarine  speed  in  knots  =  ";  V 

6570  PRINT  #1,  "distance  between  end  points  in  nautical  miles  *  ";  D 

6580  PRINT  #1,  "length  of  a  cell  side  in  nautical  miles  =  ";  SL 

6590  PRINT  #1,  "time  step  in  hours  =";  TS 

6600  PRINT  #1,  "maxiitun  number  of  time  steps  =  ";  NS 

6610  PI^INT  #1,  "track  length  in  nautical  miles  ^  ";  TL 

6620  PRINT  #1,  "delta  in  nautical  miles  =  ";  DEL 

6630  PRINT  1^^,  "leg  length  distribution  index  =  ";  NI 

6640  PRINT  #1,  "random  nuiber  generator  =  ";  R$ 

6650  IF  FLAG4  >  Q  THEN  6670 

6660  PRINT  #1,  "random  number  seed  -  ";  RNS 

6670  PRINT  #1,  "nuiber  of  tracks  completed  =  ";  NT 

6680  PRINT  #1,  "number  of  tracks  terminated  =  ";  NTR 

6690  RETURN 

6700  AS  »  NS  ♦  "  dynamic  diffusion  probability  map":  PRINT  #1,  AS 
6710  PRINT  #1,  :  PRINT  #1,  "number  of  diffusions  =  ";  ND;  PRINT  #1, 

6720  FOR  M  *  10  TO  -10  STEP  -1 

6730  XS  =  STRS(MPD(-9,  M)):  XS  =  LTRIMS(RTRIMS(XS));  PRINT  #1,  XS;  TAB(MT); 

6740  FOR  N  =  -8  TO  8 

6750  XS  =  STRS(MP0(N,  M));  XS  =  LTRIMS(RTRIMS(XS)):  PRINT  #1,  XS;  TAB((N  +  9)  *  (HT  -  1)  ♦  HT); 

6760  NEXT  N 

6770  XS  =  STRSCMPD(9,  H)):  XS  =  LTR1M$(RTRIMS(XS)):  PRINT  #1,  XS 
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;0/^W  ’NEXT  M 
iSTM  GOm  6560 
,68M(»SU8:6340- 

.6810?PRmT  #1,  "ceil  entry  sim  SM 
6820  IF  =  d;tHEMr'N  *  16  ELSE  N  =  15 
ri5830?FOR  K  «  1  TO  N 
68«b  PRiiiT  #1 , 

"6850  NEXT,  K 
6860  RETURN 

.6870  AS  »  NS  ♦  ••  static  diffusion  probability  map":  PRINT  ifi,  AS 
6880  PRINT  SI,  PRINt  #1,  "niirber  of  diffusions  =  ";  NO:  PRINT  #1, 

:6890-FOR  N  »  10  TO  -10  STEP  -1 

6900  XS  «  STRS(CIHT(P0T(-9,  N))):  XS  »  LTRINS<RTRIMS(XS»:  PRINT  #1,  X$;  TAB(HT); 

6910  FOR  N  «  -8  TO  8 

6920  XS  »  STRS(CINT(P0t(N,  M»);  XS  »  LTRIMS(RTRINS(XS»:  PRINT  #1,  XS;  TAB((N  ♦  9)  *  (HT  -  1)  +  HT> 
6930. NEXT  N 

6940  XS'»  STRS(CINT<P0T<9,  M))):  XS  «  LTRIMS(RTRIMS(XS»:  PRINT  #1,  XS 
6950  NEXT  H 
6960  GOSUB  6560 
6970,  GOSUB  6410 

'6980  PRINT  #1,  "cell  entry  sun  ■  ";  SH 
6990  PRINT  SI,  "cell  entries  are  rounded  integer" 

7000  IF  FLAG4  *  0  THEN  N  >  16  ELSE  N  «  15 
7010  FOR  X  »  1  TO  N 
7020  PRINT  SI, 

7q30.,HEXT  X 
.7040  REtURN 

7050  AS  ■  NS  ♦  "  random  tour  probability  map":  PRINT  SI,  AS 
7060  PRINT  SI,  ;  PRINT,  SI,  "nurber  of  time  steps  =  I:  PRINT  SI, 

7070  FOR  N  «  10  TO  -10  STEP  -1 

7080  XS  »  STRS(MP(-9,  M,  I)):  XI  «  LTRIM$(RTRIMt(XS)):  PRINT  SI,  XS;  TAB(MT); 

7090, FOR  N  ■  -8  TO  8 

7100  XS  «  STRKNPCN,  M,  I)):  XS  »  LTRIHS(RTRIHS<XS)):  PRINT  SI,  XS;  TAB((N  ♦  9)  *  (MT  •  1)  ♦  MT); 
7110  NEXT  N 

7.120  XS  «  STRS<MP(9,  H,  1)):  XS  «  LTRIHS(RTRIMS(XS)):  PRINT  SI,  XS 
7130  NEXT  N 
7140  GOSUB  6560 
7150  GOSUB  6270 

7160  PRINT  SI,  "cell  entry  sum  »  ";  SM 
7170  IF  FLAG4  «  0  THEN  N  «  16  ELSE  N  s  IS 
7180  FOR  X  *  1  TO  N 
7190  PRINT  SI, 

7200  NEXT  X 
7210  RETURN 

7220  REOIH  NSUMS(-LL  TO  LU,  -ML  TO  MU) 

7230  FOR  L  »  -LL  TO  LU 
7240  FOR  M  >  -ML  TO  MU 
7250  FOR  N  *  -1  TO  1 
7260  FOR  X  «  -1  TO  1 

7270  MSUMSIL,  M)  »  MTSCN,  X,  L,  M)  +  MSUMSIL,  M) 

7280  NEXT  X 
7290  NEXT  N 
7300  NEXT  M 
7310  NEXT  L 

7320  P0T(-D  /  2  /  SL,  0)  *  NT:  REM  initial  cell  probability  factor 

7330  IF  NO  «  0  THEN  RETURN 

7340  FOR  I  <  1  TO  NO 

7350  REDIM  8(-LL  TO  LU,  -ML  TO  MU) 

7360  FOR  L  »  -(LL  -  1)  TO  (LU  -  1) 

7370  FOR  M  »  -(ML  -  1)  TO  (MU  •  1) 

7380  FOR  N  *  -1  TO  1 
7390  FOR  X  *  -1  TO  1 

7400  IF  MSUMS(L  •  N,  M  -  X)  =  0  THEN  7420 

7410  B(L,  M)  =  MTS(N,  X,  L  •  N,  M  -  X)  /  MSUHS(L  -  N.  M  -  X)  *  P0T(L  -  N,  M  -  X)  ♦  8(L,  H) 

7420  NEXT  X 
7430  NEXT  N 
7440  NEXT  M 
7450  NEXT  L 

7460  FOR  L  =  -LL  TO  LU 
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'JfTO'FOR'.K  .s.^iML  TO-HU- 

'7480'roT(Li  MK*  B(L',  M);-REH  static  diffusion  probability  nap  element 

-?7500l»IEXT-X 
iTSIO’/ERASliB  -  ‘ 

>7520' NEXT  I  ■  > 

iTS'jOiERASE  MSlliiS  '  = 

'TSiJd^RETURN 

TSSd-MJS  0  " 

7560  FdR'?U- »•  -it  LU 
•7570  f.OR?Mi'«  -ML  TO  W- 
.^8o:sM  «  rotdvH)  +  SM 
75901NEXT  M 
7M0;NEXT  t 
i7610; RETURN 
7620  TrEG*«-'0 
^7630  F0R>L  '»,i  TO  NI 
7^0>  COSUB  7700- 

7650iON,ERROR;.CdTO.7660j  ;GOTO  7670 
7660  RESUME  7640 

7670  TLE6>  -DEL  *  AL  *  LOG(1  -  RAND)  /  NI  +  TLEG:  REM  trial  lefl  length  generator  for  the  relative  course  RC 

7M0,.NEXT  l,- 

769d3:RETURN 

7700:IF-RS  *  "standard"  THEN  RAND  =  RND  ELSE  GOSUB  7810:  REH  auxiliary  random  number  generator  subroutine  branch 
,  7710"=RETURN 

=7720  R$  »  "auxiliary";  REM  the  first  line  of  an  auxiliary  random  nunber  generator  setup  subroutine 
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